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We present an exploratory study of two pion photo- and electroproduction off the nu- 
cleon in the threshold region. To calculate the pertinent amplitudes, we make use of 
heavy baryon chiral perturbation theory. We show that due to finite chiral loops the 
production cross section for final states with two neutral pions is considerably enhanced. 



The experimental implications are briefly discussed. 
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I. INTRODUCTION 



Over the last few years, much interest has been focused on pion photo- and electro- 
production off nucleons. In particular, new accurate data for the processes 'jp — > n°p and 
7*p — > 7r°p close to production threshold have become available [1,2].* These have led 
to many theoretical investigations. In particular, in refs.[3] baryon chiral perturbation 
theory was used to give a model-independent description of the pertinent differential 
cross sections, multipoles and so on. While the study of these reactions, also with 
charged pions in the final state, continues on the experimental as well as on the the- 
oretical side, complementary information can be gained from the two pion production 
process 7./V — > nnN, where N denotes the nucleon and 7 the real or virtual photons. 
The two pions in the final state can both be charged, both neutral or one charged and 
one neutral. If the four-momentum of the photon is denoted by k, in case of k 2 = 
one speaks of photoproduction and for k 2 < of electroproduction. Here, we will be 
concerned with the threshold region, i.e. the photon in the initial state has just enough 
energy to produce the two pions (and the outgoing nucleon) at rest. This energy is very 
close to the first strong resonance excitation of the nucleon, the A(1232) resonance. In 
fact, presently available data focus on the resonance region and above. In that case, 
a two-step reaction mechanism of the form 7./V — > 7rA — > nnN is appropriate to de- 
scribe these data as detailed in refs.[4,5]. As we will show later, there is, however, a 
narrow window above threshold which is particularly sensitive to chiral loops, i.e. to 
the strictures of the spontaneously broken chiral symmetry. In the limit of vanishing 
current quark masses, which is a good first approximation for the two light flavors u and 
d relevant to the 7tN system, the gauge theory of the strong interactions, QCD, exhibits 
this symmetry. While QCD is formulated in terms of (confined) quark and gluon fields, 
its low-energy behaviour is dominated by the almost massless Goldstone bosons related 
to the spontaneous breakdown of the chiral symmetry. This feature allows to formulate 
an effective field theory in terms of the asymptotically observed fields (here, the pions 
and the nucleons) which is amenable to a systematic energy expansion, baryon chiral 
perturbation theory. This framework has been applied to a wide variety of hadronic and 
nuclear processes as detailed in the reviews [6,7]. In the present context, CHPT amounts 
to a calculation of tree and pion one-loop diagrams. The latter are mandated to per- 
turbatively restore unitarity and they are also demanded by the power counting scheme 
[8] underlying the effective field theory. In the baryon sector, a consistent one-to-one 

* Here, 7 and 7* stand for the real and virtual photon, respectively. In what follows, 
we will exclusively use the symbol 7 but always spell out the pertinent four-momentum 
squared which is zero for real and less than zero for virtual photons, in order. 
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mapping between the energy and loop expansion is only possible if one transforms the 
nucleon mass term into a string of vertices of increasing powers in the inverse nucleon 
mass. This can be achieved in the extreme non-relativistic limit [9]. Evidently, such a 
procedure is nothing but a series of Foldy-Wouthuysen transformations known from the 
treatment of the Dirac equation for a very heavy fermion. 

What makes the process 7./V — > 7r7riV so interesting is that, as we will demonstrate 
in detail later on, the one-loop corrections lead to a dramatic increase in the production 
rates for final states containing neutral pions. This is a counter-intuitive result since 
in the case of single pion production the cross sections for production of charged pions 
are considerably larger than the ones with neutral pions in the final state. First mea- 
surements of two-pion production at low energies have been performed at MAMI and 
we expect that the theoretical predictions discussed below will give additional motiva- 
tion to perform yet more detailed measurements of this particular reaction. Our study 
extends the one of Dahm and Drechsel [10] who discussed certain aspects of two-pion 
photoproduction in the framework of Weinberg's chiral pion-nucleon Lagrangian [11]. 

The paper is organized as follows. In section 2, we discuss the formal aspects of 
threshold two-pion production, such as kinematics, threshold multipoles and the three- 
body phase space. We also give formulae for the total cross sections for the various 
final states, in particular we write down some very handy approximations which allow 
the interested reader to get a quick and fairly accurate first estimate. We extend this 
discussion to the amplitudes close to but not at threshold since these are important in 
the chiral expansion with two charged pions in the final state (see appendix A). Section 
3 is devoted to a short discussion of the effective chiral pion-nucleon Lagrangian. We 
also discuss briefly the 'jnAN system. The latter has to be included due to the closeness 
of the A(1232) resonance. In section 4, the chiral expansion of the threshold multipole 
amplitudes is given and the various contributions are discussed. While the tree and 
one-loop diagrams in the pion-nucleon system lead to unambiguous results, the tree 
diagrams including the A have potentially large uncertainties due to certain badly known 
off-shell parameters. However, to the first two orders in the chiral expansion only one 
of these enters and the final results are found to be weakly dependent on it. Section 
5 contains the numerical results. We discuss the threshold multipoles and total cross 
sections making the approximation that the amplitude exactly at threshold represents 
the amplitude in the threshold region. For the particular channels 7^ — > tv + tv~p and 
7P — > 7T°7T°p, we extend this calculation to the first non-vanishing corrections away 
from threshold. The experimental implications are also discussed. Most of the results 
presented are for the photoproduction case since the longitudinal multipoles appearing 
in electroproduction do not carry any relevant additional information to the first two 
orders in the chiral expansion. Section 6 contains the summary and outlook. Some 
technicalities are relegated to appendix B. 
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II. TWO-PION PRODUCTION: FORMAL ASPECTS 



In this section, we will give the formalism necessary to treat two-pion photo- and 
electroproduction in the threshold region. We will only be concerned with the kinematics 
close or at threshold, the corresponding multipole decomposition and the total cross 
sections. For a more general discussion we refer the reader to [10]. 

ILL GENERAL CONSIDERATIONS 

Consider the process 7(fc) + N(pi) — > Tv a (qi) +^ b (q2) +-^(^2)7 with N denoting the 
nucleon (proton or neutron), 7 a real (k 2 = 0) or virtual (k 2 < 0) photon and pions of 
isospin a, b. The pertinent four-momenta are given (cf. fig.l) and the polarization vector 
of the photon is denoted by e M . The corresponding current transition matrix element is 

T-e= out < n a ( qi ),n b ( q2 ),N(p 2 )\r™(0)e»\N( Pl ) > in (2.1) 

with J^ m the electromagnetic current operator. From the two initial states •yp and 772 
we can form in total six final states 



IP -»■ 


7r + 7r p 


IP -»■ 


7r + 7r°r2, . 


IP -»■ 


TrVp, 


7?i — > 


7r + 7r~r2, 


7n — > 


n°n~p , 


771 — > 


yrVn. 



(2.2) 



In what follows we will mostly concentrate on the first three of these channels. To first 
order in the electromagnetic coupling e the threshold amplitudes for "fp — > 7r + 7r°?2 and 
7n — > 7r~7r°p are equal as will become clear soon. 

II.2. KINEMATICS IN THE THRESHOLD REGION 

In general, one can form five/six Mandelstam variables for the two-pion photo/ 
electroproduction process from the independent four-momenta. For our purpose, it is 
most convenient to work in the photon-nucleon center-of-mass frame. At threshold, the 
real or virtual photon has just enough energy to produce the two pions at rest. Let us 
denote by M n and m the masses of the pion and nucleon, respectively.* The threshold 
center-of-mass energy squared is 

sthr = (pi + £;) t 2 hr = (m + 2M^) 2 = m 2 (l + 4p + Ap 2 ) (2.3) 

* In section II. 5 we briefly discuss the relevance of differentiating the various masses in 
the six processes (2.2). 
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where we have introduced the small parameter 



fj, = 



m 



1 

7 



(2.4) 



The photon center-of-mass energy can be expressed in terms of s and the photon virtu- 
ally k 2 as 



ko = 



m 



2v^ 



introducing a second small parameter 



7/ 



Uthr 

K,Q — 



k 2 



1 + 2/x 



v 



fi + fi +- 



(2.5) 



(2.6) 



since we will consider only small photon four-momenta squared \k 2 \ << m 2 . When 
discussing experimental quantities like total cross sections, it is more common to consider 
the laboratory system in which N(pi) is at rest. The photon lab energy denoted by _E 7 
is then 

s-m 2 -k 2 _ (27) 



2m 



and the threshold value for two pion-photoproduction is given by 



= 2M 7T (l + /j) 



(2.8) 



For real photons only, we will calculate two-pion production amplitudes above threshold. 
In that case the pertinent four-momenta in the center-of-mass frame read 

k» = (k ,k), rf = (E 1 ,-k), p£ = (£ 2 ,-9i-? 2 ), q? = (u>i,qi), q% = M*) 

(2.9) 

together with 



q 2 + M, 



2 

7T ) 



E 2 = m 2 + kl, 



El = m 2 + q{ + ql + 2q 1 q 2 z, 



qi-k = k (u)i - xqi), q 2 ■ k = k (u 2 - yq-i, 



(2.10) 



where x,y,z denote the corresponding cosines of the angles between the three-vectors 
and ko = \k\ = (s — m 2 )/2y / i, qi = \q^\ give their lengths. From spherical geometry one 



can deduce the standard relation 

y = xz + \/ (1 — x 2 )(l — z 2 ) cos <p 



(2.11) 



with ip the azimuthal angle between the planes spanned by q± and k as well as <fi and 
q 2 . Finally, real photons satisfy e • k = 0. In the following, we will work in the Coulomb 
gauge eo = 0, so that we have in photoproduction in addition e • k = 
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11.3. THRESHOLD MULTIPOLES 

At threshold in the center-of-mass frame (i.e. <fi = q 2 = 0), the two-pion electro-pro- 
duction current matrix element can be decomposed into multipole amplitudes as follows 
if we work to first order in the electromagnetic coupling e, 

T • e = x \{io • (e x jfe) [M^ + M 2 5 ab T S + M 3 (5 aS T b + 5 bS T a )} 

+ e ■ k[N 1 5 ab + N 2 5 ab r 3 + N 3 (5 a3 T b + 5 b3 r a )] } Xl ^ 

with Xi,f two-component Pauli-spinors and isospinors and we used the gauge eo = 0. 
Clearly, for real photons only the Mi )2)3 can contribute. For virtual photons, gauge 
invariance T ■ k = allows to reconstruct T as T —T-k/ko. The multipole amplitudes 
-^1,2,3 and iV 1)2j 3 encode the information about the structure of the nucleon as probed 
in threshold two pion photo- and electroproduction. The physical channels listed in 
eq.(2.2) give rise to the following linear combination of M lj2 ,3 (and AT 12) 3 for k 2 < 0). 

7P — > 7r + 7r~p : Mi + M 2 , 

7P -> 7i + 7i°n : v^Ms, 

7P -> 7r°7r°p : Mi + M 2 + 2M 3 , 

(2.13) 

7?i — > Tr" 1 "^ n : Mi — M 2 , 

7n — > 7i°7i~p : \/2M 3 , 

7n -> 7r°7r°?i : M x - M 2 - 2M 3 

which shows the abovementioned equality for the second and fifth channel. 

11.4. TOTAL CROSS SECTIONS 

Let us consider the two-pion photoproduction close to threshold. The invariant 
matrix element squared averaged over nucleon spins and photon polarizations takes the 
form 

\M ft \ 2 = k 2 \r 11 M 1 + V2 M 2 + m M 3 \ 2 (2.14) 

with the isospin factors ?7i,2,3 given in eq.(2.13). The main dynamical assumption in 
this relation is that the two-pion photoproduction amplitude in the threshold region 
can be approximated by the amplitude at threshold. Expressing k 2 in terms of s and 
supplementing |.M/i| 2 by the photon flux factor m 2 /pi ■ k = 2m 2 /(s — m 2 ), we find for 
the unpolarized total cross section 

*l N r™ N (s) = ^{s- m 2 )T 3 (s) l^i Mi + m M 2 + m M 3 \ 2 S . (2.15) 
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Here, T 3 (s) is the integrated three-body phase space and S a Bose symmetry factor, 
S = 1/2 for the n°ir° final state and S = 1 otherwise. The integrated three-body phase 
space can be expressed as 

r 3 ( s ) = J_ [ Tl dT y / T(T + 2m)(T 1 -T)(T 2 -T) 



32tt 3 J T 3 -T 

Ti =J—(yfi- m -M nl -M 7r2 )( v ^-m + M wl + M 7r2 ), 

Ws (2.16) 
T 2 = — r (^-m-M nl + M 7r2 )(^-m + M 7rl -M 7r2 ), 

where M n i and M 7r2 stand for the masses of the final state pions and one has the 
inequality < Ti < T 2j3 . For equal pion masses an excellent approximation to eq.(2.16) 
is given by 

Of course, an analogous approximation can be derived for unequal pion masses. Conse- 
quently, the unpolarized total cross section can be approximated within a few percent 
by the handy formula 

a TN^«» {Ey) w n^l+^ lmM! + V2M 2 + m M 3 \ 2 S (E, - E^f . (2.18) 

For electroproduction, the prefactor in eq.(2.18) has to be modified slightly to account 
for the virtual photon flux normalization and then it gives the transverse total electro- 
production cross section. For the reaction — > 7r + 7r~p, we will also discuss results 
where the amplitude away from threshold is taken into account. In that case the total 
cross section has to be worked out numerically. A most efficient integration over the 
final state three-body phase space has the form 

g"tot"" 7r+7r P ( E i) = fi^ff JJ a duJ idu 2 J dx J dip\M fi \ 2 (k ,uj 1 ,u 2 ,q 1 ,q 2 ,x,y,z) 

(2.19) 

where we have chosen as the four independent variables which characterize the final state 
configuration the pion energies u>i and uj 2 , x the cosine of the angle between cfi and k 
and the azimuthal angle ip introduced in eq.(2.11). The cosine of the angle between q± 
and q 2 is already fixed by energy conservation to the value 

cw 2 - yfi(ui + u 2 ) + Ml + i(s - m 2 ) 
z = z(u 1 ,u 2 ,s) = . (2.20) 

VW- M ')K-^ 2 ) 
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Evidently, the quantity z has to lie between —1 and +1 and this condition determines 
the allowed region in the lo\U)2 energy plane. We have also considered the first above 
threshold correction for ^p — > 7t n°p as detailed in appendix A. This completes the 
necessary formalism. 

II. 5. REMARKS ON ISOSPIN-BREAKING 

In the calculations to be performed later on, it is legitimate to work with one 
nucleon and one pion mass, which we will choose to be m = m p = 938.27 MeV and 
M n = M w ± = 139.57 MeV. However, the mass difference M f ± - M v o = 4.6 MeV in 
reality leaves a 11.9 MeV gap between the production threshold of two neutral versus 
two charged pions. While we are not in position of performing a calculation including all 
possible isospin-breaking effects, a minimal procedure to account for the mass difference 
of the physical particles is to put in these by hand in the pertinent kinematics, such that 
the thresholds open indeed at the correct energy value. To be specific, for the Tc + n~p 
final state the threshold photon energy is 

E^in+n-p) = 320.66 MeV (2.21) 

whereas for tt tt°p it is 

Ef r (n°n°p) = 308.77 MeV . (2.22) 

Also, in the pertinent three-body phase space integrals we will differentiate between 
neutral and charged pion mass when we present results incorporating the correct opening 
of the thresholds. 

III. EFFECTIVE LAGRANGIAN 

In this section, we will briefly discuss the chiral effective Lagrangian underlying our 
calculation. We will also present an extension to incorporate the A(1232) resonance. 
This is mandated by the closeness of the two-pion production threshold and the location 
of the A-resonance, ttt-a — m — 2M 7r + = 14.6 MeV. Concerning the chiral interactions of 
the ttN system our discussion will be brief. Many additional details are spelled out in 
refs.[6,12]. 

III.1. CHIRAL PION-NUCLEON LAGRANGIAN 

To explore in a systematic fashion the consequences of spontaneous and explicit 
chiral symmetry breaking of QCD, we make use of baryon chiral perturbation theory (in 
the heavy mass formulation) [9] (HBCHPT). The nucleons are considered as extremely 
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heavy. This allows to decompose the nucleon Dirac spinor into "large" (H) and "small" 
(h) components 

q( x ) = e - imv - x {H(x) + h(x)} (3.1) 

with the nucleon four- velocity, v 2 = 1, and the velocity eigenfields are defined via 
yfH = H and yfh = —h. Eliminating the "small" component field h (which generates 
1/m corrections), the leading order chiral nN Lagrangian reads 

C^ N = H(ivD + g A S-u)H (3.2) 

Here the pions are collected in a SU(2) matrix-valued field U(x) 

U(x) = 4 WF 2 - tt{x) 2 + ir ■ tt{x)] (3.3) 
r 

with F the pion decay constant in the chiral limit and the so-called a-model gauge 
has been chosen which is of particular convenience for our calculations in the nucleon 
sector. In eq.(3.2) = + denotes the nucleon chiral covariant derivative, 5^ is 
a covariant generalization of the Pauli spin vector, g A — 1.26 the nucleon axial vector 
coupling constant (formally the one in the chiral limit) and 

= iu ] Vy,Uu ] (3.4) 

with u = VU and V M the covariant derivative acting on the pion fields. To leading 
order, 0(q) one has to calculate tree diagrams from 

+ ^Tv{WUV^ + X+ } (3.5) 

where the second term is the lowest order mesonic chiral effective Lagrangian, the non- 
linear cr-model coupled to external sources. At next-to-leading order 0(q 2 ) one has to 
consider tree graphs from 

C nl = H^-^D-D + ^(v-D) 2 + c 1 TTx+ + {c 2 -^)v-uv-u + csu-u 

+ C4[S", S»]u,u v - l M{ S .D,vu}- ^[S", ST] ((1 + °K v )f+ u + 5^ Tr /+ ) 

(3.6) 

Some of the terms in eq.(3.6) are the 1/m corrections from the original Dirac Lagrangian. 
But there are new terms proportional to the isoscalar and isovector anomalous magnetic 
moments k 8jV in the chiral limit or the low energy constants ci, C2, C3, C4. The latter are 
related to the tvN cr-term and ttN scattering lengths. In order to restore unitarity in 
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a perturbative fashion, one has to include (pion) loop diagrams. In HBCHPT, there 
exists a strict one-to-one correspondence between the expansion of any observable in 
small external momenta and quark masses and the expansion in the number of loops. 
In what follows we will work within the one-loop approximation corresponding to chiral 
power 0(q 3 ).* To obtain all contributions at order q 3 one has to supplement the chiral 
effective Lagrangian by additional terms C^-l + which will also serve to cancel 
the divergences of certain loop diagrams. In our case of two-pion electroproduction at 
threshold, however, all individual contributions at order q 3 are finite. 

The corrections at order g 3 to the two-pion threshold photo/electroproduction am- 
plitudes can be grouped into three classes, (i) the one-loop diagrams with insertion from 
the leading order Lagrangian eq.(3.5), (ii) kinematical corrections with inverse powers 
of m arising from the (chiral) expansion of the relativistic nucleon pole graphs and (Hi) 
contact graphs from with a priori unknown coefficients related to dynamics at 

higher mass scales and tree graphs with one vertex from C^l (e.g. the Wess-Zumino 
term which takes care about the anomalous Ward identities of QCD). In estimating the 
a priori unknown coefficients of the contact graphs we will make use of the resonance 
saturation principle and relate them to the A(1232) exchange contributions. This is 
motivated by the closeness of the A(1232) resonance and its relatively large couplings 
to the 'jnN system. The effects of heavier baryon resonances are neglected and vector 
meson exchange is suppressed by higher powers of q at threshold. Therefore we have 
now to give the chiral 7riVA Lagrangian to lowest order, 

^L = ^* a M"-( z +lh^ u H* + h - c - ( 3 - 7 ) 

with u„ = iTr(r a u v ) and we used the well known SU(4) relation g n NA = Sg^N/V^ 
together with the Goldberger-Treiman relation. The relativistic formulation using a 
Rarita-Schwinger spinor for the A introduces an off-shell parameter Z into the ttNA 
and 77riVA vertices which would be lost in the nonrelativistic isobar formulation. It 
allows us to keep track of the uncertainties in the A resonance contribution to certain 
observables. Furthermore one has to consider the 7./V — > A transition vertex which is 
proportional to the electromagnetic field strength tensor. Since the couplings of this 
vertex do not enter our final result for the two-pion production amplitudes at threshold 
calculated up to order g 3 , we do not explicitely give it here. The interested reader can 
find a discussion on this topic in ref. [14] . We have now assembled all tools to perform a 
complete calculation up to and including 0(q 3 ) of two-pion photo/electroproduction. 

* For more accurate calculations, it seems to be necessary to include all terms of order 
q 4 as demonstrated for a specific example in ref. [13]. 
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IV. CHIRAL EXPANSION OF THE THRESHOLD AMPLITUDES 



In this section, we will be concerned with the chiral expansion of the threshold 
amplitudes Mi^,3 and Ar 12 ,3. In each case we will give two complete chiral powers, the 
leading and next-to-leading term. Before giving the explicit expressions, a few general 
remarks are in order. Directly at threshold the calculation simplifies enormously. In the 
center-of-mass frame the final state nucleon and the two pions are at rest, i.e. qi = q 2 = 
0. With = (1,0,0,0) the spin-operator is the usual one = (0, a/2) and the two 
pion momenta are equal, q% = q% = (M^, 0). We therefore have in the Coulomb gauge 
eo = 0, 

e-v = e-qi = 0, S ■ qi = 0, v ■ (<?i - q 2 ) = . (4.1) 

In photoproduction, we have in addition e-k = 0. Here, the advantage of the heavy mass 
formulation (HBCHPT) clearly shows since the conditions eq.(4.1) make most diagrams 
vanishing at threshold and only very few diagrams are left to contribute. In words 
eq.(4.1) says, that the photon nucleon vertex, the (out-going) pion nucleon vertex and the 
Weinberg vertex for the two out-going pions all vanish at threshold. In photoproduction 
even the photon coupling to an out-going pion becomes zero. Furthermore, the ■yimNN 
contact term in is vanishing since it is proportional to e • v. Such selection rules 
at threshold are extremely useful to fish out the non-vanishing one-loop graphs from a 
huge number of diagrams. See also ref.[12] where such methods were applied to nucleon 
Compton scattering and ir° photoproduction at threshold demonstrating their power. 

IV. 1. THE TRANSVERSE THRESHOLD AMPLITUDES 

Let us first discuss the transverse amplitudes Mi j2 ,s- In this case we can still exploit 
the condition e • k = but we will not restrict k 2 to be zero in order to get also the 
result for electroproduction. The various contributions to Mi^.s are ordered by their 
chiral power. To leading order 0(q) we have the tree graphs from and one finds 
quite easily that due to the abovementioned selection rules all these graphs are exactly 
zero at threshold. The leading nonzero contributions therefore comes from tree graphs 

(2) 

with one insertion from C^^. The two non-vanishing diagrams are shown in fig. 2. It 
is interesting to note that the first non-vanishing term for two-pion photoproduction 
at threshold start at order q 2 in all channels. This is very different from single pion 
photoproduction where the charged channels are strongly enhanced. 

Next, we have to study all contributions at chiral power 0(q 3 ). Among these are the 
one-loop contributions. Again the selection rules eq.(4.1) and e-k = reduce drastically 
the number of graphs which are nonzero at threshold to four. These are depicted in 
fig. 3. Finally we have to work out all polynomial contributions at order q 3 . Part of 
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these (the 1/m 2 corrections) can be most easily obtained if we calculate the nucleon 
pole graphs relativistically and then expand in powers of M n and k 2 . Further terms 
are given by tree graphs involving two vertices from C^ N where at least one carries a 
constant c±, C2, C3, C4 or k s ,k v . The two posssible graphs are shown in fig. 4a, their sum, 
however, vanishes since the nucleon propagators in both have opposite sign. Then we 
have to study tree graphs with one vertex from C^l. Most of these insertions just lead 
to a multiplicative renormalization of the 0(q) graphs, which are, however, all equal 
to zero. C^l contains the Wess-Zumino term incorporating the anomalous (natural 
parity violating) vertex 7 — > 3n. Nevertheless, the diagram shown in fig. 4b vanishes at 
threshold, simply because it is proportional to e^apq^ = with qi = q2- Finally, in 
order to complete the list of all possible contributions at 0(q 3 ) we have to consider the 
(genuine) contact terms from C^^. These carry new low energy constants which are not 
known a priori and have to be determined from phenomenology. We will invoke here the 
resonance saturation principle (which has been tested in the meson sector to work very 
well) and estimate these constant from the A(1232)-resonance contribution to two-pion 
photo/electroproduction at threshold. One expects sizeable effects from the A(1232) 
since first it is quite close to threshold and second its couplings to the 771 N system are 
very large (about twice the nucleon couplings). On first sight the distance of only 14.6 
MeV of the A(1232) from threshold seems to give rise to overwhelming contributions 
since one naively expects that the very small denominator 



m A ~~ s thr — m — 2-Mtt ttia + m + 

enters the result. We have worked out all possible graphs with a single or double A 
excitation (cf. fig. 5) and found that at threshold the small and dangerous denominator 
tiiA — m — 2M n always gets cancelled by exactly the same term in the numerator. This is 
a very important feature since otherwise the M n expansion would be in serious trouble. 
Mtt-/ (mA — m — 2M 7T ) formally counts as small, 0(q), but numerically it is of course very 
large. After this detailed discussion of the chiral expansion of the transverse threshold 
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amplitude, let us now give the final result, 



M 1= ^ + <% 2 ), 



4 m 2_p2 



Mo 



4mF 2 

2 



(2^ -!-«„) + 



4m 2 F 2 



, eg A M n / 8 + 4r —- r 
H — r < ; arctan v 1 + r 



+ i 



y/3(2 + r) 1 + r + r z 



1 + r 



(l + r )3/2 



In 



1 + r + r 2 

" (l +r )3/2 

2 + r+ v/3(l + r) 



B A 



7t r 
h arctan - 

2 V / TT7 



M 3 =^^(l + ^-2^ 2 ) + 



8mF 2 



8m 2 F 2 



+ 



\/l + r + r 2 

eg 2 A Mir 



lemm^F 2 



2r 



+ — 777 o , arctan V 1 + r 



1 + r + r z 



+ i 



VS(2 + r) 1 + r + r 2 2 + r + ^3(1 + r) 



1 + r (1 + r) 3 / 2 
with the ratio r = — k 2 /AM 2 and 



Vl + r + r 2 



1 + r (1 + r) 3 / 2 
+ 0(q 2 ) 



7T r 
h arctan . 

2 v^T7 



= 2mi + mAm " m2 + 4Z[m A (l + 2Z) + m(l + Z)] 
m A — w 



(4.3) 



(4.4) 



which involves the off-shell parameter Z of the 7riVA vertex. In fact, taking the allowed 
range of Z given in ref.[14], i.e. — 0.8 < Z < 0.3 we find a weak Z-dependence, i.e. 9.9 
GeV < £? A < 15.1 GeV. Furthermore, from the isospin factors of eq.(2.13) we see that 
to order M n the A contributions are absent in the 7r°7r° channels. 

It is interesting to note that Mi is zero to O(q ) and 0(q). At next order q 2 it 
only receives a contribution from the relativistic iV-pole graphs, but no loop or counter 
term contribution. The result eq.(4.3) has been written completely in terms of phys- 
ical parameters, not the chiral limit values which enter the effective Lagrangian. For 
F^.M^.gA and m these differences only would show up at the next order 0(q 4 ). The 
isovector anomalous magnetic moment, however, has been renormalized and includes a 
non-analytic piece ~ Vt% [12], 



9A mM ^ 

4ttF 2 



(4.5) 



Note that the first terms in the expansion of arising at order q 2 differ from 
the corresponding expressions given in ref.[10] by the term proportional to k v . In that 
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reference, the Weinberg Lagrangian was used and the photon was coupled in via minimal 
substitution. Such a procedure can, however, not generate the photon nucleon coupling 
proportional to the anomalous magnetic moment. It is furthermore important to write 
down such an anomalous photon nucleon vertex in a manifestly chiral invariant fashion 
using in eq.(3.6) the quantity and not just the photon field strength tensor. When 
expanding /+ in powers of the pion field the pertinent 'ynnNN vertex (proportional to 
k v ) based on chiral symmetry is automatically generated. 

Another point worth mentioning is that the transverse multipoles are /in- 

dependent only through their loop contribution. This can be understood from the fact 
that the tree graphs have to be polynomial in both M w and k 2 and that a term linear 
in k 2 is already of higher order in the chiral expansion. It is also not possible to further 
expand the r-dependent functions since r = — k 2 /AM 2 counts as order one and all terms 
have to be kept. 

At first sight the multipoles M 2) 3 seem to behave singular in the chiral limit = 
since then r becomes infinite. This is, however, not the case and the chiral limit is 
perfectly smooth (as it should be) with 

O loop o loop eCl\ I — e rr, OH, 

M 2 =-2M 3 = I J^V^=^[^-F 2 (fc 2 )] (4-6) 

ov 

where F 2 [k ) is the nucleon isovector magnetic form factor in the chiral limit calculated 
up to order q 3 (see ref.[12]). 

The loop contribution to the transverse multipoles of two pion production as given 
in eq.(4.3) have a nonzero imaginary part even at threshold. This comes from the 
rescattering type graphs. Due to unitarity the pertinent loop functions have a right hand 
cut starting at s = (m + M^) 2 (the single pion production threshold) and these functions 
are here to be evaluated at s = (m + 2M 7T ) 2 (the two-pion production threshold). For 
photoproduction the complicated fc 2 -dependence of the loop terms disappears and we 
have the complex constants 

M i° P = ff|f { | + *[2v/3 - ln(2 + VS)} } 



(4.7) 



This completes the discussion of the transverse threshold multipoles up-to-and-including 
order q 3 . 
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IV.3. LONGITUDINAL MULTIPOLES AT THRESHOLD 



In the electroproduction case, we also have the longitudinal threshold amplitudes 
iV"i 5 2,3. Since we can no more exploit the condition e • k = the photon coupling to 
an out-going pion line is non- vanishing and therefore we obtain a nonzero contribution 
already at leading order O(q) involving a pion propagator (cf. fig. 6a). Adding up all 
terms which arise at order q and q 2 we find the following results 

N x = 0(g) , 

N2 ~ F 2 (4M 2 ~k 2 ) + AmF 2 + ° {q) ' (4-8) 
N 3 = ~N 2 + 0(q). 

It is interesting to note that none of the low energy constants Ci, 02,03,04 and the 

(2) 

anomalous magnetic moments which enter C show up in the final result. The graph 
of fig. 6b contains the higher derivative tttiNN vertex proportional to ci, 02,03 but the 
isospin factor 5 ac e c3b + (a <-> b) = annihilates this diagram. With eq.(4.8) we have given 
two full chiral powers (leading and next-to-leading) for the expansion of the longitudinal 
threshold amplitudes AT 1;2 ,3. We do not work out here the third chiral power, which 
comes from loops and polynomial counter terms. We will now return to photoproduction 
and present our results for observables like total cross sections. 

V. RESULTS AND DISCUSSION 

In this section, we will first present numerical results for the threshold multipoles 
and total cross sections in the isospin limit. For the experimental implications, we also 
consider the kinematical shifts for the various channels as discussed in section 2.5. 

V.l. RESULTS IN THE ISOSPIN LIMIT 

First, we must fix parameters. We work with m = m p , M n = M n ± : = 93 MeV, 
g 2 N /ATx = 14.28 and e 2 /4-7r = 1/137.036. In fig. 7, we show the invariant matrix-elements 
\Mij I for the 7p initial state versus the off-shell parameter Z defined in eq.(3.7) (k 2 = 0). 
As stressed after eq.(4.4), the A does not contribute to the 7r°7r° final state. For the 
allowed range of Z taken from ref. [14] , the Z-dependence is weak. The surprising result 
is the dominance of the 7r°7r° final state for the threshold multipoles chirally expanded 
to 0(M 7r ). In table 1, we give the 0(1) and 0(1) + 0(M 7r ) contributions with and 
without the A term. To lowest order, the 7r°7r° final states are suppressed, as already 
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noted in ref.[10]. It is interesting to compare these results to the case of single pion 
photoproduction. In that case, the Kroll-Ruderman term dominates charged production 
and the loop effects are small. For -fp — > jv°p, however, the loop corrections are quite 
large [3]. The difference to the two pion production case is that the total cross sections 
for photoproducing a charged pion are nevertheless much larger than the corresponding 
ones for neutral pions (in the threshold region). Notice also that for the n + n~ final 
state the 0{M 1T ) loop and kinematical corrections are largely cancelled by the ones from 
the A. In case of electroproduction, the k 2 -dependence of the Mj (i = 1, 2, 3) is weak. 
For k 2 ranging from to -0.1 GeV 2 , the \M.ij\ (in GeV -3 ) change from 11.2 to 9.4, 
20.1 to 15.1 and 35.2 to 29.7 for -fp — > tv + tv~p, 7r + 7r°n and 7r°7r°p, respectively. The 
corresponding longitudinal multipoles iVj (i = 1,2,3) are completely dominated by the 
pion pole (compare eq.(4.8)) and are not shown here. 

In fig. 8, we show the total cross sections for the first three channels given in (2.2) 
by approximating the amplitudes in the threshold region through their threshold values. 
For E 1 = 330 MeV, we find* 

0tot(7P 
0tot(7P - 

These are very small cross sections, but we refer to section 5.2. for further discussion. 
In fig. 9, we compare the lowest order cross section for the production of two charged 
pions with the first correction from O(q), cf. appendix A. While the former goes like 
(E 1 - E^) 2 , the latter approximatively rises as (E 1 - E^f. At E 1 ~ 323 MeV, these 
two contributions are of the same size. This narrows the window in which the 7r°7r° 
final state is dominant to the first few MeV above threshold. The first correction to 
this particular channel is rather small as discussed below (cf. appendix A) . Let us now, 
however, review these results in the light of the different thresholds already mentioned 
a few times. 

V.2. EXPERIMENTAL IMPLICATIONS 

To connect to the experimental situation, we now consider the three-body phase 
space with the physical masses for the corresponding pions. This automatically takes 
care of the various threshold energies. In the loops we work, however, with one pion 

* These numbers differ from the ones reported in ref.[15] since in that reference the 
A-contribution was omitted. 



> tv tv°p) = 0.36 nb, 

Tr+TrV) = 0.22 nb, (5.1) 
n + n-p) = 0.08 nb . 
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mass since that effect is small as discussed in appendix B. Chieffy, the loop functions 
are sensitive to (M n + — M 7T o)/M 7T + = 0.03, i.e. a few percent effect. Fig. 10 shows the 
calculations with the correct phase-space and using the threshold matrix-elements. For 
7P — > 7t + tt~p, we also show the first correction above threshold from At E 1 = 320 

MeV, the total cross section for 7r°7r° production is 0.5 nb whereas the competing 7T°Tv + n 
final state has utot = 0.07 nb. Double neutral pion production reaches 0"t o t = 1.0 nb 
at = 324.3 MeV in comparison to a tot (7P — > 7r°7r + ?i) = 0.26 nb and a to t(7P —> 
n + n~p) < 0.1 nb. This means that for the first 10. ..12 MeV above 7r°7r° threshold, one 
has a fairly clean signal and much more neutrals than expected. Remember that to 
leading order in the chiral expansion, the production of two neutral pions is extremely 
suppressed. Of course, the above threshold correction for this channel, which comes from 

(2) 

£)rN ( an d higher orders) should be calculated systematically. The first correction, which 
vanishes proportional to |^| (i = 1,2) at threshold, has been calculated (see appendix 
A) and found to be very small. The corresponding cross section at E 1 = 320, 325 and 
330 MeV is a t fi r f corr = 0.009, 0.026 and 0.056 nb, i.e a few percent of the leading order 
result. It is, therefore, conceivable that the qualitative features described above will not 
change if even higher order corrections are taken into account. 

VI. BRIEF SUMMARY AND OUTLOOK 

In this paper, we have performed an exploratory study of two-pion photo- and 
electroproduction off the nucleon. We have focused on the threshold region and used 
heavy baryon chiral perturbation theory to calculate the first two powers in the chiral 
expansion of the transversal (Mi^) and longitudinal (Ar 1)2 ,3) mulitpoles defined in 
eq.(2.12). The pertinent results (we focus on the photoproduction case) are: 

(i) To leading order in the chiral expansion, only the multipoles M 2 and M3 are non- 
vanishing, with M 2 = — 2M 3 . Therefore, the production of two neutral pions is 
strictly suppressed. 

(ii) At next order in the chiral expansion, one has to consider kinematical (1/m) correc- 
tions, one-loop contributions and tree graphs with insertion of the A(1232) (chirally 
expanded). The loop contributions make the 7r°7r° final state dominant in contrast 
to the expectation from single pion production. 

(iii) We have calculated total cross sections for the various channels eq.(2.2) under the 
assumption that the amplitudes in the threshold region can be approximated by 
the exact threshold amplitude. Taking into account the different pion masses in the 
respective phase spaces, these cross sections are shown in fig. 10. There is a window 
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of about 10 MeV above 7r°7r° threshold in which one can detect much more neutrals 
than expected. 

(iv) We have considered the first above threshold corrections for the processes ^p — > 
n + n~p and 'jp — > n n°p. While the former are sizeable even very close to threshold, 
the latter are rather small. We conclude that ultimately one should consider the 
full amplitude above threshold (including also more dynamics from the resonance 
region) to draw more quantitative conclusions. We hope to report on this in the 
not too far future. 



APPENDIX A: LEADING ORDER AMPLITUDE ABOVE THRESHOLD 

In section IV, we considered the two-pion photoproduction amplitude exactly at 
threshold and we have seen that the selection rules make the 0(q) contribution exactly 
vanishing at threshold in all channels. Here we will now work out the amplitude for 
the process 'jp — > n + n~p (with k 2 = 0) away from threshold to leading order O(q). 
For the 7r°7r° channels to this order the amplitude above threshold is still equal to zero 
everywhere . The five diagrams shown in fig. 11 (plus their crossed partners with the 
pion lines interchanged) give rise to the following amplitude, with the kinematics taken 
over from eq.(2.9), 



T ■ e = — - < -e-q\ - e • q 2 r H a ■ q 2 a ■ e H a ■ e a ■ q\ 

2F4 y qi • k q 2 -k u) 2 ui 

0J 2 qi ■ k uiiq 2 ■ k J 



(Al) 



To calculate the unpolarized total cross section we have to average over initial spins and 
sum over final spins, i.e. we must compute the quantity 

^Tr[T • 6 (T • c)t] (A2) 

where the trace is taken over 2x2 spin-matrices. From (A. 2) we get a quadratic form 
in the polarization vector e which we have to average over the two photon polarizations 
perpendicular to k. The pertinent polarization averages are done via the substitutions 

(f • 9l) 2 - i<h 2 (l - A, 

(r-sa'-^d-A (A3) 

e-qie-q 2 -> -qiq 2 (z-xy), 
(?) 2 -l. 
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with x,y,z the cosines of the inclined angles. Therefore, the invariant matrix element 
which enters the cross section integral eq.(2.19) takes the form 

I kA |2 e2 / 9l^2 /-, 2x , o ^152^1^2 , v . ql^l n 2 , 

o 2 [ 2/1 2x ^2 ( 1 , q2(k y - qiz)\ t 2/i 2 ^ wi /M , qi(k x-q 2 z) 
~ l 9a Qi I 1 — x ) 

+ qMz-xy)(-^-(±+ Ql{k0X _-? Z) ) q2(k ° y " QlZ) 
+ 9a 



<?i ■ ft \o;i u 2 qi ■ k J q 2 ■ k \uj 2 uj\q 2 ■ k 



qi ■ k \u 2 uiq 2 ■ k J q 2 ■ k \Ui u 2 qi ■ k 

2«? 24 qlqtjl-x*) 2 _ , ggg 2 2 (l - y 2 ) 2 _ 

W 2 + a; 2 + u/ 2 ( ? i-A0 2 ( 01 + 9lJ + u/ 2 (<6 • ife) 2 ( 



H 2<?lg2 : (&o(> + yz - 2xy 2 ) + q 2 (xy - z)) H 2glg2 . (fc (y + - 2x 2 y) + ^(xy - , 

<^1^2<?2 • « ^1^2<?1 ' « 

4^192 . 2qlql(z-xy) , 2 Jl 

xy H f fc (2xy - z) - /c <?iy - fe <?2^ + <?i<?2j f 

a>ia> 2 ^i^2<?i • kq 2 - k K J J 

(AA) 

The lengthy expression of the invariant matrix element from the rather simple amplitude 
(A.l) makes clear that calculations above threshold are rather voluminous. In fact, 
for the first 10 MeV above threshold the last two terms in eq.(A.l) could be safely 
neglected since they vanish with a higher power of \qi\ (i = 1, 2). In order to represent 
the most general amplitude above threshold more than a dozen of amplitude functions 
are needed to characterize a specific channel (see ref.[10]) and each of them depends 
on five independent variables (e.g. E^, u>i, oj 2 , x, y). In comparison to this enormous 
complexity single pion photoproduction is rather simple. The process is completely 
described in terms of four amplitude functions which depend on E 1 and e.g. the cms 
scattering angle. 

Let us briefly turn to the first above threshold correction for — > 7r°7r°p. The 
corresponding tree diagrams must include at least one insertion from C^ N as noted 
before. From these, the ones with the -firNN vertex from C^ N go to zero as \q t \ (i = 1, 2) 

(2) 

as one approaches threshold whereas the ones with a photon-nucleon vertex from 
and two pion-nucleon vertices from C^ N go as \q x \ \q 2 \, i-e. they vanish faster. So to 
leading order in small momenta, we consider only the first type of diagrams and find 



T.e= 1691 



2mFl 



— cr • (e x q 2 ) H a ■ (e x q x ) 

L0 2 Uli 



(A5) 



and from this the total cross section follows as 

o2„4 r r /, ,2^2 , ,2^2 



a - W = 384.3^4 J J z2<i **i **» + -5- + 2q iq2 zj . (A6) 
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In this case, one can perform the x and (p integrations in |.M/i| 2 analytically. We have, 
however, also performed the four-dimensional integration numerically as a check and 
found good agreement with the results based on (A. 6). 

APPENDIX B: LOOP FUNCTIONS 

The calculation of the one-loop graphs shown in fig. 3 leads to the following loop 
integrals (in d space time dimensions): 

7 / (27c) d v ■ l(M 2 - l 2 )(M 2 _ (/ + fc )2) = ^73(w, k 2 ) + 

7 / {2-kY v ■ l(M 2 - /2)(M2 - (/ + k) 2 )(M 2 -(l + k- W) = 9 ^ Qs{uJ ' M "' ^ + ' ' ' 

(B.l) 

with uj = v ■ k and q = M T v. All propagators are understood to have an additional in- 
finitesimal negative imaginary part in the denominator, which makes the Wick-rotation 
to euclidean space time unique. The ellipsis stand for terms proportional to v^v u , k^k v 
etc., which are not needed. Using standard Feynman parameter representations we 
obtain for the symmetric sums at d = 4 which enter the threshold amplitudes 

1 r 1 

^(u,k 2 ) +-f 3 (-u,k 2 ) = — / dx^M%-x 2 u 2 +x{x-l)k 2 , 
Qsiu.M^.k^+Qsi-u.-M^.k 2 ) 



x 



327rM w (w - M x ) 

f dx{^Ml - x 2 u 2 + x(x - l)k 2 - y/M 2 + AxM^u - M n ) - x 2 u 2 + x(x - l)k 2 }. 
Jo 

(B.2) 

These are then to be evaluated at uj = 2Af 7r + zO. Let us finally comment on isospin 
breaking effects. If we consider 7r°7r° production the pions in the loop are all charged, 
therefore M n = M n + in eq.(B.2) and the functions are to be evaluated at u = 2M 7r o +i0. 
Since both real and imaginary part of the functions in eq.(B.2) are smooth around 
uj = 2M 7T + with finite first derivative the isospin breaking effect from the different pion 
masses is of order (M n + — M^o)/M^+ ~ 3% and therefore very small. This is quite 
different to the reaction -fp — > iv p where the respective loop functions vary strongly 
around the branch point uj = (of square root type) and the effective isospin breaking 
parameter is v / (M 7r + - M 7r o)/M 7r + ~ 18%. 
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FIGURE CAPTIONS 



Fig.l The process ^(k) + N(p 1 ) — > n a (qi) + 7r 6 ((/ 2 ) +N(p 2 ). Here, 7, n and N denote the 
photon, the pion and the nucleon, in order, a, b are isospin indices, e is the polar- 
ization vector of the photon and the corresponding four-momenta are indicated. 

Fig. 2 Lowest order diagrams which lead to the leading order result for the multipoles 

(2) 

Mi 2 3. The circle —cross denotes an insertion from jC nN . 

Fig. 3 One-loop diagrams contributing to M 1)2 ,3 at order M n . 

(2) 

Fig. 4 (a) Tree diagrams with two insertions from C j^. (b) Tree diagram sensitive to the 
anomalous 7 — > 37r vertex. These vanish as described in the text. 

Fig. 5 Classification of the tree diagrams with A(1232) insertions, (a) and (b) lead to the 
lowest order 0(M n ) contribution, (c) e.g. contains the dangerous energy denomi- 
nator as described in the text. Crossed diagrams are not shown. 

Fig. 6 (a) Lowest order diagrams contributing to the longitudinal multipoles AT 1)2 ,3. (b) 
The circle-cross denotes an insertion from C^^. The terms proportional to q 
(i = 1, 2, 3, 4) and k v>8 do not contribute as described in the text. 

Fig. 7 Invariant matrix-elements (in GeV -3 ) for the ^p initial state (k 2 = 0) as a function 
of the off-shell parameter Z (in the isospin limit). 

Fig. 8 Total cross sections (in nb) for the ^p initial state (k 2 = 0) as a function of the 
photon lab energy E 1 (in the isospin limit). 

Fig. 9 Total cross sections (in nb) for the 'jp — > n + n~p. The dashed line gives the result 
from the exact threshold amplitude and the solid one the first above threshold 
correction (see appendix A). 

Fig. 10 Total cross sections (in nb) for the ^p initial state (k 2 = 0) with the correct three- 
body phase space. The various thresholds are indicated. 

Fig. 11 Tree diagrams with insertions exclusively from C^ N which lead to the first above 
threshold correction for -fp — > tv + tv~p, eq.(A.l). 
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TABLE 



process 


0(1) 


0(l) + 0(M 7r ) 


0(l) + 0(M n ) 








without A 


with A 


IP - 


■> TV + TV~p 


1 1 ni 




1199 


IP - 


-> TV + TV°n 


7.78 


14.18 


20.06 


IP ~ 


-> 7r°7r°p 


0.00 


35.25 


35.25 


771 - 


■> 7T + 7r _ 72, 


11.01 


12.84 


10.72 


771 - 


-> 7r~7r°p 


7.78 


14.18 


20.06 


771 - 


-> 7r°7r°n 


0.00 


30.96 


30.96 



Table 1: Invariant matrix-elements \r\\M\ + 772M2 + 773M3I for the various two 
pion production channels in GeV -3 (k 2 = 0). 
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